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Universal crossovers and critical dynamics of quantum phase transitions: 
A renormalization group study of the pseudogap Kondo problem 
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The pseudogap Kondo problem, describing a magnetic impurity embedded in an electronic envi- 
ronment with a power-law density of states, displays continuous quantum phase transitions between 
free and screened moment phases. In this paper we employ renormalization group techniques to 
analytically calculate universal crossover functions, associated to these transitions, for various ob- 
servables. Quantitative agreement with the results of Numerical Renormalization Group (NRG) 
simulations is obtained for temperature-dependent static and zero-temperature dynamic quantities, 
at and away from criticality. In the notoriously difflcult realm of finite-temperature low-frequency 
dynamics, usually inaccessible to both NRG and perturbative methods, we show that progress can 
be made by a suitable renormalization procedure in the framework of the Callan-Symanzik equa- 
tions. Our general strategy can be extended to other zero-temperature phase transitions, both in 
quantum impurity models and bulk systems. 

PACS numbers: 75.20.Hr,74.70.-b 



I. INTRODUCTION 

Continuous quantum phase transitions^ are a subject 
of intense research, both in experimental and theoretical 
condensed matter physics. Particularly interesting are 
the low-temperature, low-energy properties near quan- 
tum critical points, which are crucially influenced by 
the non-trivial properties of the quantum critical ground 
state. A theoretical description requires the knowledge 
of the low-energy long- wavelength effective theory for the 
critical degrees of freedom. Critical exponents, charac- 
terizing the power-law behavior of observables upon ap- 
proaching the critical point, may then be obtained via 
suitable renormalization group (RG) techniques. A re- 
lated, but considerably harder, task is the calculation of 
full crossover functions describing the universal behav- 
ior of observables, e.g., along the changes that take place 
from the high-temperature quantum critical region to the 
low-temperature regimes, which appear in the vicinity 
of simpler (ordered or disordered) quantum mechanical 
ground states. These crossover functions are particularly 
relevant for experiments, as they allow to analyze mea- 
surements in terms of data collapse and universal scaling. 

On theoretical side, a large body of work has con- 
centrated on quantum critical points above their upper 
critical dimension where observables can be calculated 
in straightforward perturbation theory, but where uni- 
versality is limited due to ultraviolet (UV) divergencies 
and the absence of hyperscaling. In contrast, interacting 
quantum critical points below their upper critical dimen- 
sion require a more sophisticated treatment, and rela- 
tively little is known regarding full crossover functions. 
Exceptions are models for insulating magnets, with many 
results summarized in Ref. 1. 

An interesting class of model systems showing quan- 
tum phase transitions, with both interacting and non- 
interacting critical points, are quantum impurity mod- 



els. The transitions here are conceptually simpler than 
in bulk systems, as only a finite number of degrees of free- 
dom becomes critical at the transition point. ^ Powerful 
numerical techniques, most notably Wilson's Numerical 
Renormalization Group'^ (NRG) , have been developed to 
study the asymptotic low-energy regime of impurity mod- 
els, and thus allow to benchmark analytical approaches. 
Besides being model systems for fascinating many-body 
physics, including the Kondo effect and local non-Fermi 
liquid behavior,** quantum impurity models are of great 
relevance for nanofabricated quantum dots as well as for 
correlated bulk quantum systems treated within dynam- 
ical mean-field theory*^ (DMFT) and its generalizations. 

The so-called pseudogap Kondo problem^'^ displays 
particularly rich physics including various impurity quan- 
tum phase transition. Here, a quantum spin 1/2 couples 
to conduction electrons with a semi-metallic density of 
states (DOS) that vanishes as p{uj) ex \ujY {r > 0) around 
the Fermi level. Although only special values of the pseu- 
dogap exponent r can be experimentally accessed to date 
(namely r = in metals, r = 1 in graphene layers and d- 
wave superconductors, and r = 2 in p-wave superconduc- 
tors), this model in its generality displays a challenging 
variety of non-trivial properties. Recently, it was real- 
ized that all critical fixed points of the pseudogap model 
can be accessed perturbatively after having identified the 
values of r that correspond to (upper or lower) criti- 
cal dimensions of the problem. ^'^ Subsequently, critical 
properties of the model have been analytically calculated 
using perturbative RG techniques. The results for vari- 
ous critical exponents and static observables at criticality 
(such as the entropy or the anomalous Curie constant) 
show impressive agreement with NRG simulations.^ This 
approach illuminates the need to consider the adequate 
degrees of freedom to describe a given quantum phase 
transition (see Ref. 10 for a simple illustration). 

The purpose of this paper is to extend the analytical 



approach of Ref. 9 to the calculation of full crossover 
functions for static and dynamic quantities, both at 
and away from criticality. Technically, the calcula- 
tion of the scaling behavior of the physical quantities 
will be achieved by using the powerful Callan-Symanzik 
equations^^"^'^ within the field-theoretic RG, allowing us 
to compute the universal crossovers functions of the most 
interesting physical observables within the renormalized 
theory. We will restrict our attention to the particle- 
hole symmetric case and to the range of the DOS ex- 
ponent r = ... 1/2, but our techniques can be applied 
to the particle-hole asymmetric situation^ near r = 1 
as well. Comparison with numerical data from NRG 
will be made, with extremely good agreement regarding 
both the temperature dependence of static observables, 
and the zero-temperature dynamic correlation functions. 
With respect to the finite-temperature dynamics, it is 
known that direct perturbation theory displays unphysi- 
cal singularities in the low- frequency regime.^'* We show, 
however, that the use of Callan-Symanzik equations may 
allow some progress in obtaining physically sound re- 
sults. Because the regime where frequency is smaller 
than temperature is also notoriously hard to access via 
NRG simulations,^^ our results could possibly be useful 
in the future for benchmarking the quality of this nu- 
merical approach. In a more general context, we believe 
that the powerful methodology employed here (being rel- 
atively unfamiliar in the condensed matter context) will 
also be of interest for general studies of quantum phase 
transitions, e.g., in quantum impurity models or bulk 
magnets. 

We note that some scaling functions to be determined 
below can also be obtained within the large- A'^ expansion 
of multichannel Kondo models, both for the metallic^®'^^ 
and pseudogap^* density of states. In addition, a partial 
study of the dynamics at the critical point using renor- 
malization ideas was undertaken in Ref. 19, and the local 
moment approach^'' was also applied to this problem and 
compared to NRG.^^ We will comment on these results 
in the body of the paper. 

The bulk of the paper is organized as follows: In Sec. II 
we briefly review the physics of the pseudogap Kondo 
model, in particular the possible quantum phase tran- 
sitions and the associated observables. Focusing first on 
the weak-coupling regime of the pseudogap Kondo model, 
which is perturbatively controlled for small values of r, 
we develop in Sec. Ill the necessary formalism to compute 
the full universal crossovers functions of the observables, 
and compare the analytic results to NRG data. In Sec. IV 
we address the regime where r is not small, using a previ- 
ously developed effective theory,^ valid when r is close to 
the value 1/2. The calculation of the universal crossovers, 
employing the language of the Anderson model, parallels 
the procedure used with the Kondo model, and is again 
compared with success to the numerical data. A field- 
theoretic derivation T-matrix in the Kondo problem is 
given in the appendix. In the course of the paper, some 
previously unaddressed technical issues will also be dis- 



cussed, such as the absence of double-counting problems 
in renormalized perturbation theory. 



II. THE PSEUDOGAP KONDO PROBLEM: 
SUMMARY OF PREVIOUS RESULTS 

A. Kondo and Anderson models 

The pseudogap Kondo problem^ (for a review, see 
Ref. 2) originates from the question of how a mag- 
netic impurity behaves when coupled to a fermionic bath 
with a depleted (power-law) single-particle spectrum at 
the Fermi level. This question is of interest to various 
condensed-matter systems such as d-wave superconduc- 
tors, graphene, and zero-gap semiconductors (and is in 
contrast to s-wave superconductors, where there is a hard 
gap). The qualitative understanding of the pseudogap 
Kondo problem is that Kondo screening, i.e., the quench- 
ing of the local degrees of freedom, is weakened with re- 
spect to the metallic case due to the lack of low-energy 
states in the bath. This results in the possibility to have 
quantum phase transitions controlled by the strength of 
the Kondo coupling, the amount of particle-hole symme- 
try breaking and the shape of the DOS. In the following, 
the DOS will be fixed at the simplified model form 
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which has an algebraic dependence in energy character- 
ized by the number r > (the standard metallic Kondo 
problem corresponds to r = 0). Here, e^ are single- 
particle energies of the fermionic bath, D denotes the cor- 
responding half-bandwidth, and Nq = {r + l)/(2i?''+^). 
We consider the Kondo model with an antiferromagnetic 
coupling J > between the impurity and the bath: 



Hk = E^'»4.c,. + '^^-E4(o)^v(o) 



+ £;o5]4(0)c,(0) 



(2) 



in standard notation. For _Eo = the model is particle- 
hole symmetric, and Eq tunes the degree of asymmetry. 
A related model is the pseudogap generalization of the 
Anderson Hamiltonian, which reads: 

Ha ^ J2'^''^l'^^ka + Ud\d^dld^ 

k<7 

+ Y.i^ddUa + VldicM + H.c]). (3) 

a 

It is worthwhile noting that the Anderson model is equiv- 
alent to the Kondo model in the limit of frozen charge 
fluctuations, i.e., single occupancy. This so-called Kondo 
limit of the Anderson model, is obtained via Schrieffer- 
Wolff transformation"* in the limit of large U , and leads to 



the effective parameters of the Kondo Hamiltonian (2): 
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The particle-hole symmetric case corresponds thus to the 
value Erf — —U/2. 



B. Fixed-point structure 

One of the main results of the extensive numerical 
studies performed by Gonzalez-Buxton and Ingersent^ 
was the generic phase diagram of the pseudogap Kondo 
model (as well as of the pseudogap Anderson model, be- 
longing to the same universality class). A subset of the 
results, for an intermediate value of the bath exponent 
r, is displayed in Fig. la as a function of Kondo cou- 
pling J and particle- hole asymmetry Eq. At particle-hole 
symmetry, Eq = 0, a, quantum critical point (denoted 
by SCR) separates the local-moment (LM) phase from 
the strong-couphng (SSC) regime. The SCR fixed point 
is expelled from the decoupled LM point as soon as r 
is non zero, but collapses with the SSC point when r 
reaches the value 1/2. Away from particle-hole symme- 
try, the transition is controlled by SCR only in the range 
< r < r* ~ 0.375, and for larger values of r, a sec- 
ond critical point (denoted ACR because of its intrinsic 
asymmetric nature) drives the transition. Fig. la thus 
corresponds to the intermediate regime r* < r < 1/2 
where both critical points coexist. Recently, a general 
understanding of this quite complex phase diagram has 
been reached on an analytical level, ^ thanks to the iden- 
tification of the effective theories which can be used in 
the various regimes of the parameter r. 

For small values of r, the Kondo model is the correct 
effective model to describe the phase transition between 
the local-moment and the strong-coupling phases. In- 
deed, the critical SCR point at J — Jc moves towards 
the LM fixed point as r goes to zero, and the quan- 
tum phase transition at particle-hole symmetry can be 
controlled by (renormalized) perturbation theory in J 
within the Kondo Hamiltonian (2), as discussed in pre- 
vious works. ^'^ It is important to point out that this ex- 
pansion corresponds to the situation of a lower critical 
dimension, where the critical fixed point is located adja- 
cent to the trivial fixed point around which the expansion 
is carried out (here the "ordered" LM phase), in analogy 
to the renornialization flow of the non-linear sigma model 
above two dimensions. ^'^ In principle, we would like to 
compute any physical quantity along the crossover from 
SCR to LM using the RG technique - this will be done 
in Sec. III. 

On the other hand, when r approaches the value 1/2, 
the critical SCR point moves to strong coupling and the 
Kondo Hamiltonian is no longer amenable to a pertur- 
bative treatment. In this case, an expansion around the 




SSC 



FIG. 1: (color online) Generic flow diagrams for a) the Kondo 
model (with r* < r < 1/2); b) the particle- hole symmetric 
Anderson model (with r < 1/2). In both cases, the universal 
crossovers (depicted as colored lines) will be computed. See 
the main text for a description of the nomenclature used in 
this figure. 



SSC fixed point is more natural^, and the particle- hole 
symmetric Anderson model (3) provides the adequate ef- 
fective theory, with perturbation theory in U now being 
controlled in the small expansion parameter e = 1 — 2r. 
Again, the fiow diagram in Fig. lb refiects a situation of 
lower critical dimension. This full crossover will be cal- 
culated in Sec. IV. It is important to note that - within a 
perturbative framework - we cannot describe crossovers 
associated to a runaway flow of the coupling constant. 
Therefore, the crossovers SCR — > SSC for small r and 
SCR — > LM for small (1/2 — r) small are inaccessible to 
analytical calculations. 

In the particle-hole asymmetric case, not to be ad- 
dressed in this paper, the critical fixed point ACR for 
r > r* can be accessed using the perturbative approach 
of Ref. 8, whereas for r < r* multicritical behavior ob- 
tains, controlled by the SCR fixed point. 



C. Observables 

To establish some notation, we will introduce the ob- 
servables that we will consider in this work, and summa- 
rize some results previously obtained in Ref. 9. 



1. Impurity susceptibility 

Magnetic susceptibilities are generically obtained by 
coupling an external field to both the bulk electronic de- 
grees of freedom and the impurity part: 



Afield = / dxK{x) 



\ci{x)^c^,{x) + hi^p ■ S (6) 



where S is the impurity spin [it corresponds to 

X^crcr' ^a'^^^f^^a' ™ ^^^ '"^^^ ^^ ^^"^ Anderson model]. The 
bulk field h-^ varies slowly as function of the space co- 
ordinate, and /limp is the magnetic field at the location 
of the impurity. With these definitions, a spatially uni- 
form field applied to the whole system corresponds to 
hu = /limp = h. Response functions can be defined 
from second derivatives of the thermodynamic potential, 
fi = — TlnZ, with respect to both hu and /limp- Thus 
Xu,u measures the bulk response to a field applied to the 
bulk only, Ximp,imp is the impurity response function to 
a field which is applied to the impurity only, and Xu.imp 
is the crossed response of the bulk to an impurity field. 

Following standard practice, we define the impurity 
contribution to the total susceptibility (associated to a 
spatially uniform magnetic field applied both to the bulk 
and to the impurity) as 



Ximpi-i j — X 



imp, imp 



I ^Xu,imp + (,Xu,u Xu,u Jj \') 



where x^^lt is the susceptibility of the bulk system in 
absence of the impurity. 

For an impurity spin of size S = 1/2 we expect the 
spin to remain unscreened in the whole LM phase: 



Ximp(r^ 0, J < Jc 



1 



(8) 



(Here and in the following we will use units such that 
n = kB^ 1.) 

At criticality Ximp does not acquire an anomalous di- 
mension (in contrast to Xiocj see below), because it is the 
response function associated with a conserved quantity, 
the total spin of the system Stot- In fact, one finds a 
renormalized Curie law. 



Ximp(7'^ 0, J= Jc) = 



(^imp('') 
f 



(9) 



which can be interpreted as the Curie response of a frac- 
tional spin, since the universal number Cimp(i') < 1/4. 

Finally, the strong-coupling fixed point SSC provides 
a modified (albeit trivial) Curie response: 



Xi 



,p(r^O,J> J,) = — . 



(10) 



All these results are summarized in Fig. 2. 

The phase transition occurring at particle-hole sym- 
metry is described by an interacting fixed point and thus 
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FIG. 2: (color online) Low-temperature limit of TXimp (con- 
tinuous upper line) and Txioc (dashed lower line) as a function 
of the Kondo coupling J. 



obeys hyperscaling. The static susceptibility must there- 
fore satisfy the relation 



1 ( T 

Ximp(,J J "^ j = 7p;*Pimp I 7^ 



T 



T* 



(11) 



where T* ex \J ~ JcY" is the energy scale above which 
quantum critical behavior is observed, and v is the as- 
sociated correlation "length" exponent. In this formula 
^imp is a universal crossover function. Such a scaling 
form (and the other ones written below) will be derived 
in Sec. III. 



2. Local susceptibility 



The local susceptibility is given by 



Xloc(r) =X: 



imp, imp : 



(12) 



which is equivalent to the zero-frequency impurity spin 
autocorrelation function. In the LM phase, one has a 
Curie law xioc oc 1/T at low temperature, with a prefac- 
tor not pinned to the value 1/4, as, contrarily to Ximp, 
the local susceptibility is not associated to a conserved 
quantity. We can then define a "magnetization" rriimp as: 



lim Tx\oc{T,J < Jc) 

T-.0 



imp 



(J) 



(13) 



which turns out to be a suitable order parameter for the 
quantum phase transition,^ as this quantity vanishes con- 
tinuously when J -^ Jc, see Fig. 2. Physically, the mo- 
ment of size 1/2 is spatially smeared, and mimp measures 
its fraction which remains localized at the impurity site. 



At criticality, m 



imp 



0, and the local susceptibility 



follows a non-trivial power law. 



X\oc(T ^ 0, J = Jc 



T'^-Vx 



(14) 



defining the universal anomalous exponent 77^ which con- 
trols the anomalous decay of the two-point correlations 
of the impurity spin (see below) . One obtains a similar 
behavior for xioc at J > Jc, although 77^ takes then a 
trivial value. 

Because of the hyperscaling property, the static local 
susceptibility follows 



xioc(r, J) 



Bi 



rpl 
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(15) 



where $1 is a universal crossover function, whereas Bi is 
a non-universal prefactor (this is related to the fact that 
Xioc is associated to a non-conserved quantity). 

Hyperscaling also implies a scaling in uj/T for the dy- 
namic quantities. The dynamic local susceptibility can 
be defined from the imaginary-time correlation function: 
X\oc{t,T, J) = {S{t) ■ 5'(0)). In particular, the spectral 
density of its Fourier transform obeys the following scal- 
ing form 



Xioci^^T,J) = 



B2 



jl-'J; 
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uj T 



(16) 



with again $2 a universal crossover function, and B2 a 
non-universal number. 



3. Impurity entropy 



In the remainder of the paper, we intend to perturba- 
tively calculate all of the scaling functions defined above. 



III. CROSSOVER FUNCTIONS NEAR r = 0: 

WEAK-COUPLING REGIME OF THE KONDO 

MODEL 

As discussed in Ref. 9, the Kondo model is the correct 
starting point to describe the phase transition between 
the LM and SSC phases when the parameter r is small. 
In analogy to the theory of critical phenomena, r = con- 
stitutes a lower critical dimension of the problem, and the 
RG for the pseudogap Kondo model is analogous to the 
RG in the non-linear cr-model. The RG equations have 
first been derived by WithofF and Fradkin,^ who pointed 
out the possibility to have a quantum phase transition 
driven by the strength of the Kondo coupling. 

One can switch to the dimensionless renormalized cou- 
pling j = NoZ^ Zf^^J, where /i is an arbitrary renor- 
malization energy scale, and Zj (Zj) are the vertex (field) 
renormalization factors. To one-loop order, the RG equa- 
tion reads6'22.26.27 



m 



dln/i 



rj - J 



(19) 



and the renormalization factors evaluate to 



Quantum critical points in impurity models show in 
general a finite residual entropy (this is to be contrasted 
with bulk quantum critical points, where the entropy has 
to vanish at zero temperature). The flow of the coupling 
constant in the vicinity of the critical point will then im- 
ply a crossover of the entropy, in a similar fashion to the 
behavior of the Curie constant associated to the impurity 
susceptibility. 



4- T-matrix 

An important quantity in quantum impurity models 
is the conduction electron T-matrix, which describes the 
scattering of conduction electrons off the impurity. For 
the Anderson model, the T-matrix is given by the relation 
T{uj) = V^Gd{i^), with Gd{uj) the full d-level propagator. 
For the Kondo model, one defines a propagator Gt of the 
composite operator To- — S ■ ^^, c^, (O)Tcrcr' such that the 
T-matrix is given by T(cj) = ,PGt{'-^) (see App. A for a 
derivation of this formula, as well as Refs. 22-25). 

The hyperscaling ansatz gives the general expression 



r(c.,r,J) = Src.'-''-$r(^,|^ 



(17) 



where an anomalous exponent 777- was introduced. One 
has in fact the exact result^^ rjr — 2r, so that the critical 
T-matrix behaves as 

r(u;,r = 0,J= Je)cx^. (18) 



z, = ^ + -. Zf = l 



(20) 



The beta function (19) yields an infrared unstable fixed 
point at 



3c = r + 0(r2) 



(21) 



which controls the transition between LM and SSC. Be- 
cause j is dimensionless, we will prefer to write the dif- 
ferential equation (19) w.r.t. a dimensionless parameter 
A, and integrate it between the values j for A = 1 (corre- 
sponding to the initial value of the renormalized coupling 
at the renormalization scale /x) and j{X) at the lower ob- 
servation energy scale A/i, with A < 1. The solution 
reads: 



m = 



l + [{r-j)l]]\-^ 
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(22) 
(23) 



We have introduced the energy scale T* , which van- 
ishes at the quantum critical point, and generally sets 
the crossover scale between the various fixed points. We 
notice first that the running coupling obeys on the one 
hand j(A) ^ r ii j ~ r (fixed-point behavior), while 
j(A) oc A'' if j < r and \^ <C T* (flow to the weak- 
coupling LM fixed point). In this case, the running cou- 
pling constant is therefore always perturbatively small 



+ < 



FIG. 3: Diagrammatic expansion for the impurity suscep- 
tibility of the Kondo model, with the tree level contribu- 
tion (of Ximp.imp type) plus the order j perturbative correc- 
tion (of Xu.imp type). Here full lines denote the Abrikosov 
fermion propagator, while dashed lines represent the bulk 
fermions. Each vertex, labelled by a dot, yields the value 



(when r ^ 1), and one should be able to follow the full 
crossover between the critical point and the LM phase. 
On the other hand, for j > r the above expression shows 
a pole, that signals that the crossover between the criti- 
cal point and the SSC phase occurs through runaway flow 
from weak coupling, invalidating the perturbative expan- 
sion. This regime is analogous to the metallic Kondo 
problem, where the perturbative expansion breaks down 
at the scale Tk, the Kondo temperature,^ which can be 
here associated to T* for j > r (one can check that as 
r vanishes, the usual one-loop expression for the Kondo 
temperature is recovered from T*). In the following, we 
will focus on the perturbatively accessible crossover be- 
tween SCR and LM (corresponding to j < r), while the 
crossover from SCR to SSC will be tackled through the 
mapping onto the Anderson model discussed in Sec. IV. 



A. Regular observables 

Within this section we will calculate the crossover func- 
tions for observables which do not show anomalous power 
laws. Technically, this means that the renormalization Z- 
factors for these observables are strictly unity and pertur- 
bative corrections do not change the power-law behavior 
of the tree level result, but imply a renormalized univer- 
sal amplitude. Usually such observables are associated 
with conserved quantities, ^^ such as the total spin of the 
system in the case of the impurity susceptibility. 



1. Impurity susceptibility 

One can calculate the crossover function for the im- 
purity susceptibility from the perturbative result derived 
in Ref. 22. The diagrammatics is done introducing a 
constrained Abrikosov fermion /J to represent the im- 
purity spin S — X^ctct' /1^^/ct" '^it^ *^^ constraint 
^^/T/o- = 1. To lowest order in j, the impurity sus- 
ceptibility is given in Fig. 3. The correction of order j 
takes the form 
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which, when evaluated in the limit of infinite UV cutoff 
D and at lowest order in r, gives 
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(25) 



using the tree-level result ximp = 1/{4:T) at j ~ 0. An 
intuitive way to obtain the full crossover from SCR to 
LM is to set /i = T in the previous expression and to 
plug in the running coupling j(A), Eq. (22), taken at a 
scale A = T/ ji. However, we would like to derive this 
result rigorously in the simplest case of Ximp, in view of 
the more complicated physical quantities we will consider 
later on. 

In the following we resort to the concepts of 
renormalization, ^^'^^ familiar in the ficld-thcorctic con- 
text. We start by noting that the bare impurity suscep- 
tibility xSip explicitly depends on the high-energy cutoff 
D but should be insensitive to the arbitrary energy scale 
/i, while conversely the renormalized Ximp is independent 
of D (which is ultimately taken to infinity in the renor- 
malized theory). The two are in general related through 



xtpiT,J,D) 



•^impXimp(J j .?j Mj 



(26) 



where the renormalization factor Zjmp = 1 here, because 
one is dealing with a conserved quantity free of multi- 
plicative renormalizations. We therefore have 



l^-^xEnpiT,j,D) 
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(27) 
Ximp(T,j>) = (28) 



with (3{j) = dj/dln/i. This is a special case of a Callan- 
Symanzik (CS) equation, which reflects the invariance of 
the renormalized theory with respect to the high-energy 
cutoff, and therefore is a strong statement of universality. 
It is solved by: 



Ximp(T,j,^) =F In^ 



/3(J') 



(29) 



where F is an arbitrary function. From the /3 func- 
tion (19), one has InA = J dj' / j3{j'), which inserted 
in the previous equation gives: 

Ximp(r,j,M) = FhniXfi)-!' ^j (30) 

=^ Ximp(T,j,^) = Ximp(T,j(A),A^) (31) 

We finally use the fact that T, fj,, and 1/ximp have di- 
mensions of energy, to obtain: 



1 



Ximp(T,j,^) = -Ximp ( -,j(A),^ 



T 



(32) 



This type of equation is the starting point for the cal- 
culation of the crossover functions defined in Sec. II C. 



Because fi appears identically on both sides of (32), we 
will drop it in the following. We proceed by inserting 
the perturbative result (25) in the r.h.s. of the previous 
equation: 



ap[J^ T J) — ,rp 



T 
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m+o{f 



(33) 



We now discuss the choice of the renormalization pa- 
rameter A which is crucial for all subsequent computa- 
tions. Generally, an equation like (33) is valid for any A 
(provided that all terms of the expansion are included!), 
however, for the application of a perturbative scheme we 
have to require the higher-order corrections to be small. 
As these in general have a singular structure, there is 
usually a unique choice of A (up to a perturbative re- 
definition) to obtain a controlled perturbative expression. 
In the present case of Ximp, P corrections (displayed ex- 
plicitly in Sec. IIIB below) display ln(T/A/x) singular- 
ities. This fixes A to A = T / ^, to allow perturbation 
theory to be valid, and this yields: 



Ximp(T,j) = ^[l-j(T^)] 



(34) 



Using the running coupling j(A), Eq. (22), we finally get: 



J^ Ximp 
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1- 
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(35) 



which is the desired scaling form (11). Formula (35) in- 
terpolates smoothly between the free-spin value 1/4 at 
low temperature {T <^ T*) and the non-trivial fractional 
spin value (1/4) (l — r) associated to the quantum critical 
high-temperature regime (T 3> T*). 

Expression (35) can be compared to NRG simulations 
(see Figs. 4 and 5), with T* as the unique fitting param- 
eter. For r = 0.1 the agreement is seen to be excellent, 
over 40 orders of magnitude! We note that for the larger 
value r = 0.3, the fit deteriorates, but only at high tem- 
perature. This is because the critical point does not lie 
anymore in the perturbative regime (a two-loop calcu- 
lation may improve this result). At low temperature, 
however, the agreement becomes better, a fact that we 
interpret from the flow of the running coupling constant 
towards weak coupling. 
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^imp 




FIG. 4: (color online) Impurity susceptibility rximp(r) at 
r = 0.1 given in Eq. (35) (solid), capturing the crossover 
between SCR at high temperature and LM at low temper- 
ature. Also shown are NRG data (dots) for Tximp, ob- 
tained for a pseudogap Anderson modeP" with U/D = 1, 
NqV^/D — 0.05. The crossover scale T* has been fitted: 
T* /D ~ 1.7 X lO"'^'^. [The reason for the very slow crossover 
from SCR to LM is the large value of the correlation length 
exponent, v — 1/r + 0(1) ~ 10.] Note also that deviations 
from scaling behavior occur when temperature is comparable 
to the bandwidth D, as seen by a kink in the last data point. 
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2. Impurity entropy 

The same calculation in the case of the impurity en- 
tropy shows the perturbative result :^^ 



FIG. 5: (color online) Impurity susceptibility rximp(r), 
Eq. (35), now for r = 0.3, together with NRG data for 
U/D = 1, NoV^/D = 0.25, resulting in T* = 2 x 10"^^ 
As discussed in the text, the agreement between the analyt- 
ical curve and the NRG data is better at low temperature, 
where the running j(A) has flown towards weaker coupling. 
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Application of the CS equation with A = T/ ji gives: 
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(37) 



Since the observed variations in the impurity entropy are 
actually quite small for the pseudogap Kondo model, ^ we 
will not try to compare this quantity with the numerical 
results. 



B. Local susceptibility 



In this section we will present a detailed calculation of 
the static and dynamic local susceptibility, xioo within 
the weak-coupling expansion of the Kondo model us- 
ing the CS equations. In contrast to the quantities 
of Sec. Ill A, xioc does display singular behavior with 
anomalous power laws. 

The starting point is the calculation of the renormal- 
ization factor Z^, since we now deal with a quantity 
which is not conserved. The corresponding diagrammatic 
expressions are in Refs. 26,27, and yield [this result will 
be re-derived below, see Eq. (44)]: 



The factor /^ (A) is at the origin of all anomalous power- 
law behaviors as we will see shortly, and reads: 



/xW 



AC^'"^' 




(43) 



At the quantum critical point, j = r, the above expres- 
sion reduces to /^(A) = A'' . 

Note that the choice A = lo/ ^ in (42) proves rigorously 
the uj/T scaling property stated in Eq. (16). It is im- 
portant to understand that such scaling behavior gener- 
ally applies only below the upper critical dimension, but 
as the particle-hole symmetric Kondo problem only pos- 
sesses two lower critical dimensions, the scaling ansatz is 
fulfilled for all values of < r < 1/2 (outside this range 
the critical point ceases to exist). In the particle- hole 
asymmetric case, hyperscaling is violated when r > 1 due 
to the presence of a dangerously irrelevant coupling.^ In 
such a case, the observable vanishes in the naive scal- 
ing limit (because the fixed point value of the coupling 
is zero), and the leading behavior is exposed only upon 
considering corrections to scaling (which do not follow 
hyperscaling laws). 



Z^ = l 



3_ 
2r 



Oif 



(38) 



As before, one can relate the bare susceptibility X\lc *° 
the renormalized one xioc"^^ by: 



Xi^cC'^; T, J, D) = Z^x\oc{i^. T,j, n) 



(39) 



The next step is to derive the bare quantity with respect 
to the renormalization scale fi, which after some straight- 
forward algebra leads to the Callan-Symanzik equation 
for the renormalized local susceptibility 



,^|+^(-^')|+^-(^') 



Xioc(^,T,j,m)=0 (40) 



with the anomalous exponent 



f + Oi,f). 



(41) 



[Because Zi^p = 1, the CS equation (28) for Ximp lacks 
such an extra term.] The main relation in the following 
calculations is the integrated Callan-Symanzik equation 
for the local susceptibility, which is given by 



Xioc(c^,T,j) = i/^(A)xioc(ppj(A)) (42) 



1. Impurity magnetization 

A quantity, which turns out to be a good order pa- 
rameter for the transition is the impurity magnetiza- 
tion, which is related to the low-temperature limit of 
the zero-frequency local susceptibility by Eq. (13). We 
start, as previously, by a direct perturbative calculation 
of xioc(r,j) = xioc(w = 0,T,j). To second order in J 
the diagrams shown in Fig. 6 evaluate to: 
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dx I dy\xy\y\ 
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2- - 1 - e^'-y 
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Note that the quantity in brackets is the perturbative ex- 
pression for the hare susceptibility, where the bare cou- 
pling J has been expressed in terms of the renormalized j, 
and the limit of infinite bandwidth has been taken. Fur- 
ther, the additional renormalization factors Zj and Zf 
do not appear as they can be approximated by unity to 
the order we are working here. Evaluating the integrals 
yields: 
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FIG. 6: Diagrams contributing to the local susceptibility xioc 
of the Kondo model up to order j^, with conventions given 
in Fig. 3. For the second and third graphs above, the self- 
energy-like contributions require a mass substraction. 



where 7 — 0.5772... is Euler's number. 

We can now apply the CS equation (42) at zero fre- 
quency: 



Xioc{T,j) = jfxWx\oc[T^JW 



(45) 
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FIG. 7: (color online) Rescaled static local susceptibility 
TxiociT) for r — 0.3 according to formula (46). Although 
not performed here, a comparison to NRG data could be pos- 
sible. 



We see again that the particular choice A = T/fx al- 
lows to control the perturbative expansion. Putting to- 
gether (43) and (45), we therefore arrive to the expres- 
sion: 



Xioc{T,j) = Y 






exp 



l + \T/T*\ 



with Ci 



r - J 



1 — 1" In —— 



(46) 



which is the scaling form (15), with the anomalous expo- 
nent 



Vx 



(47) 



In the limit T ^ T*, we find, as expected, that the 
magnetization vanishes as a power law at the quantum 
critical point: 



Tx\oc{T,j) = m^ 



imp 



with mimp oc {jc - jY^^ 



(48) 
(49) 



to perform the perturbative calculation on the imaginary- 
time axis. The graphs displayed in Fig. 6 give a com- 
plicated integral (not shown), which, however, simplifies 
considerably by taking a double time derivative: 
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(51) 
(52) 



where (3 above defines the inverse temperature 1/T. This 
expression is readily integrated, using Eq. (44) to obtain 
the integration constant, which gives the simple result: 

X1oc{t) == 1-J In . (53) 

TTl 

Now performing the Fourier transform and the analyt- 
ical continuation, one finds: 



However, when T ^ T* , critical behavior with an anoma- 
lous temperature dependence obtains: 



Xioc{T,j) ex 



1 



Ti- 



(50) 



as anticipated in Eq. (14). The complete crossover func- 
tion, defined in equation (46), is shown in figure 7. 



2. Dynamic susceptibility 

We are now interested in the imaginary part of the real 
frequency susceptibility Xioc('^i^j j)- It is actually easier 



Xloc('^) 
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e^'fi 
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Trf3ujS{uj) + vr— . (54) 
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The zero- frequency contribution corresponds to the Curie 
part of the static susceptibility, see Eq. (44), while the 
finite-frequency term is associated to dissipation. Let 
us first focus on the zero-temperature case. Applying 
the CS equation (42), one gets at finite frequency (for 
readability, we neglect here the exponential factor which 
only provides a small rcnormalization of the amplitude) : 






Xfi 


2r 
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(55) 
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ation in the limit w ^ T, as was already noted in the 
context of quantum antiferromagnets by Sachdev.^^ 

We will now show that the application of the CS 
equation before the analytic continuation (namely on 
the imaginary time axis) does allow to recover a phys- 
ically sound result (although the general applicability 
of this procedure remains to be benchmarked) . In- 
deed, the perturbative expression (53) implies the choice 
A — 7rT/[/isin(7rTT)] of the renormalization parameter 
to remove the logarithmic divergency. For simplicity, we 
assume that one sits at the critical point [j — r), so that 

fxW = ^^ i'^ equation (43). The CS equation for the r- 
dependent susceptibility gives immediately (up to a small 
multiplicative correction) : 



FIG. 8: (color online) Rescaled dynamic local susceptibility 
uJx'i'oci^jT = 0) at zero temperature for r = 0.1, with the 
same parameters as in Fig. 4. NRG data points are the dots, 
and the solid line follows formula (56). Since the crossover 
scale T* was already determined from Ximp(r), the only fit- 
ting parameter here is the non-universal amplitude C2. 



Employing A = w//i (this is again vindicated by requir- 
ing that the next term appearing in perturbation theory 
should be smaller than the present leading contribution) , 
we finally obtain: 



UJ 



1 + \uj/T*\- 



T* 



(56) 
This scaling function crosses over from ui ^+'^^ behavior 

at w <C T* to w~ "'"'' behavior at uj ^ T*, as can be seen 
in Fig. 8. 

We have thus obtained the zero-temperature limit of 
the scaling function (16), and can now wonder whether 
this program can be fulfilled at finite temperature as well. 
Of course, since the regime where to ^ T obeys the zero- 
temperature scaling form, one is rather interested in the 
so-called "relaxational" [to <^ T) and crossover (w ^ T) 
regimes. It is easy to see that a direct application of the 
renormalization equations on the real frequency axis fails 
for u) <^T . Indeed, by the Kramers-Kronig relation, the 
static local susceptibility is related to the spectral density 
by: 



X\oc{T) = duj 



XIA^,T) 



(57) 



At non-zero temperature, the l.h.s. is finite and given by 
Eq. (46), so that x't'oci^'"^)' ^^ *^^^ function of lu, has 
to vanish at low frequency. However, the 1/lu divergency 
shown in (54) is in contradiction with the above state- 
ment [note that the application of the CS equation cannot 
allow to fully remove this l/w divergency]. The reason for 
the failure of this calculation is the non-commutativity of 
the renormalization procedure and the analytic continu- 
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(58) 



This well-defined function can be Fourier-transformed 
and analytically continued:^^'^^ 



Xioc(w,T,j = r) = - 



t\ ii ) ^ 



- sin(7rrV2)r(l - r^) 



F(rV2 -I- iuj/2TTT) 
F(l - r2/2 - iuj/2TiT) 



(59) 



One can then expand this expression with respect to lo 
and j = r, which gives two different results according to 
the order in which the limits are taken: 



lim lim xioc(w,r,j = r) = - + i — — ;j, 

j2 

hm hni xioc(w, T,j = r) ^ n— . 
w— >0j— >o UJ 



(60) 



The first expression is physically correct, but is clearly 
non-perturbative, while the second expression simply re- 
covers the singular perturbative correction. 

This explicit calculation has clearly demonstrated the 
non-commutativity of the epsilon-expansion and the ana- 
lytic continuation, but has hinted that progresses can be 
made when one proceeds in the correct order, namely by 
the use of the CS equation on the imaginary-time pertur- 
bative expression, followed by the analytic continuation. 
In principle, the computation of Xioc('^i^jJ) away from 
the critical point can be done along this line, although a 
numerical analytic continuation would be required. 



C. T- matrix 

We will now turn to the calculation of the scaling func- 
tion for the T-matrix of the pseudogap Kondo problem, 
along the crossover from the quantum critical into the 
local-moment regime. Besides the comparison to NRG 
data, an important question here is again the determi- 
nation of the dynamical behavior, including the low- 
frequency relaxational regime. We will see that some 
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control of the perturbation theory can be achieved due to 
the pecuUar logarithmic nature of the finite-temperature 
corrections in the frequency-dependent T-matrix. 

The T-matrix is identically zero at the decoupled LM 
fixed point, and starts at second order in the Kondo cou- 
pling, with the result: 



(61) 



T(°)(zt^„) = I sgnK)|c^„r J^iVo / da:- 



Note that the bare coupling J appears here. This implies 
that this prefactor will not follow an RG fiow within the 
renormalized CS formulation we are using, and simply 
acts as an non-universal prefactor. For simplicity, we 
will omit the prefactor J^A'o in the intermediate steps 
of the following calculation, it will be re-instated in the 
final result. We now want to calculate the full crossover 
function for this quantity, and start by setting up the 
Callan-Symanzik formalism. 



+./ 




FIG. 9: Diagrammatic expansion of the T-matrix of the 
Kondo model up to next-to-leading order (J'^), with conven- 
tions given in Fig. 3. 



1. Callan-Symanzik equations forT 

The basic starting point is as previously^ ^ 
T''iLo,T,J,D) = ZTTiLo,T,j,^l) 



(62) 

where T^ denotes the bare T-matrix, and T the renor- 
malized T-matrix. The renormalization factor is^^ 



Zt = l-2j /r + 0{f) (63) 

at one-loop order. This gives 

VrU) = m^^^ = -2j + 0{f) . (64) 



A general derivation of Zj- is also given in App. A, which 



shows the relation Zr = Z- ^Z?, so that the bare cou- 



pling can also be expressed as J = /i^''^.^^ ' Nq^j. Since 
fidJ/dfi = 0, one gets the relation rjrij) — —2r + (3{j)/j. 
Although rjxij) picks up contributions to all orders in 
perturbation theory, the existence of a critical point jc 
such that /3(jc) = implies^^ the exact result 777- = 
VtUc) = — 2r. This leads to T{u!) ~ l^l^^^ behavior 
at the quantum critical point, as we will see below. 

The diagrammatic expression for the T-matrix at order 
J is shown in Fig. 9 including all spin indices. On a 
technical level it is crucial to perform the perturbative 
expansion for the propagator of the complete operator 
Ta = E.'S-T^.'C,,, e.g. T^ = (1/2)(/|A - /[/, )c. + 



ft 



/ /| C| . If only the last term of the expression is chosen, 

namely fJf^c^, one is actually computing a correlation 
function which is not associated to a physical observable, 
and one encounters double counting problems, which lead 
to an incorrect result for Zr- 

The same steps already followed to obtain the CS equa- 
tions for the local susceptibility give now: 



r(c.,T,j) = /^(A)A'-T(^,^,j(A) 



(65) 



using that the bare scaling dimension of T is r (hence 
the A'' prefactor). The anomalous behavior is provided 
by the function: 



/r(A) = cxp 



"^dA' 



TVr ijiX')) 



{T*y + \Xfi\\ 



(66) 
(67) 



which reduces to /r(A) = A^^*" for j = r. We now fol- 
low the RG philosophy by calculating the r.h.s. of equa- 
tion (65) in renormalized perturbation theory. 

The diagrammatic expansion seen in Fig. 9 gives the 
result for the imaginary part of the T-matrix at order j: 



-T"iu;,T,j) = \u;\^ 



l-2jJdx\x\^P-^ 



exp(-a;|w|/r) + 1 



eix] 



(68) 



where P denotes the principal value. In this formula, the 
pole coming from the Zt factor was taken into account, 
and appears as the 9 term inside the integral. For w = /i, 
it regularizes the expression in the UV so that Tpert (^ = 
/i) is finite. However, to make sense of this formula when 
Lu ^ ^, one has to do a small r expansion (up to order 
r"), leading to a cancellation of all l/r poles, and to the 
result: 
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exp(-a;|w|/r) + 1 



- e{x) (70) 
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(for consistency, we have also expanded the \x\'^ term to 
lowest order). 

Using the CS equations, we will investigate in turn 
the zcro-tcniperature crossover regime (T = 0, j < r), 
and then the finite-temperature critical regime (T > 0, 
j = r), the latter both at large and small frequency, 
corresponding to critical and relaxational dynamics, re- 
spectively. 



2. Zero-temperature scaling function for T 

At zero temperature, the integral in Eq. (70) vanishes, 
and one is left with a single frequency-dependent loga- 
rithmic divergency, similar to the logarithmic temper- 
ature singularity in the perturbative expansion of the 
static local susceptibility, formula (44). The applica- 
tion of the CS equation (65) clearly imposes the choice 
A = a;//x to eliminate this large term, which, using 
Eqs. (65-67), gives: 



-r"(w,r = 0,j<r) = J^iVo 



W'>2r 



[(T*)- 



(71) 



where we have re-introduced the non-universal prefac- 
tors. This function crosses over from uj'^ / {T*)'^'^ behavior 
at w <C T* to tt'"'' behavior at w ^ T*, as can be seen 
on Fig. 10, and displays an excellent agreement with the 
NRG data. 



T"(a;, T = 0) 




FIG. 10: (color online) T-matrix -T"{uj,T = 0) at zero tem- 
perature for the same parameters as in Fig. 4. NRG data 
point are the dots, and the solid line follows formula (71), 
where the non universal amplitude J'^NofJ.'^'' has been fitted 
to match the numerics. Note that the high-frequency peak at 
uj ~ U/2 in the numerical result corresponds to charge fluc- 
tuations (Hubbard bands) of the Anderson model employed 
in the NRG simulations, and cannot be accounted for in our 
calculation based on the Kondo model. The inset shows the 
same curve on a double-logarithmic scale. 



3. Finite-temperature scaling function for T 

For simplicity, we first focus on the finite-temperature 
behavior of the T-matrix above the quantum critical 
point, j = r. At T = 0, or when lu '^ T, the previous 
arguments apply and the choice of the renormalization 
scale A — uj/ij, together with the CS equation (65) leads 
to the critical behavior: 



-r'{iu,T,j 



1 



r oc 



for w > T 



(72) 



Note that the integral in (70) is of order 1 for uj >T and 
provides only a small temperature correction to the lead- 
ing critical behavior, even when frequency is comparable 
to temperature. 

We now investigate the low-frequency behavior of the 
finite-temperature T-matrix. As discussed before (see 
Ref. 14 and Sec. IIIB2), one generally obtains large 
terms in the expansion from perturbation theory, pre- 
venting a controlled calculation. 

Examining the integral (70) for tj <C T, one clearly sees 
that a ln(ci;/T) divergency is pulled out, and combines 
with the already present logarithm to give: 



-T"(c^,T,j) = |c.r[l + 2jln(T/^) + 



(73) 



where . . . denote regular to /T corrections at small fre- 
quency. This shows that indeed large terms are generated 



in the perturbative expansion in the relaxational regime. 
However, the extra singularity is logarithmic, which sug- 
gests that the RG could be used to resum it. Indeed, 
applying the GS equation (65), using the above pertur- 
bative result on the r.h.s., and taking into account that 
/r(A) = A"^'' at the quantum critical point, we get: 



-T"{uj,T,j = r) 



A2' 



[l + 2rln(r/A/x) 



(74) 



Now clearly the good choice for the scale A is 

A = T/fi, (75) 

so that perturbation theory is controlled. This in turn 
implies the following small uj behavior: 

-r\u;, T, J = r) « ^ for ^ « T (76) 

up to small perturbative corrections. We have to re- 
mark here that, although the co^ frequency behavior for 
T{uj,T,j = r) at Lo <^ T is reminiscent of the zero- 
temperature result near the local moment fixed point, 
i.e., Eq. (71) for T{lj,T = 0,j < r) at cj < T*, criti- 
cality is reflected in the anomalous temperature depen- 
dent prefactor in Eq. (76). However, the fact that both 
frequency dependencies are identical appears to be the 
reason behind the success in computing the relaxational 
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FIG. 11: (color online) T-matrix -T"{oj,T,j) for r = 0.3 
at finite temperature from the formula (79) obtained via the 
analytic RG. Different curves correspond to three values of 
j < r, associated from top to bottom to T* /T — 0, 1, 100 
(the topmost curve corresponds thus to the quantum critical 
point). Note that no NRG data in the regime uj < T is 
available to date. 



of the finite-temperature scaling function is very rem- 
iniscent of the result obtained from the local moment 
approach, ^'^ and a quantitative comparison to these re- 
sults could be interesting. We note that our analyti- 
cal estimate Eq. (76) should apply to the multichannel 
pseudogap Kondo model at large- A^ as well.^^ As the T- 
matrix does not saturate for oj < T, but rather shows 
the behavior (76), we conclude that the conformal scal- 
ing ansatz introduced for the metallic case^^ does not 
hold in general for the T-matrix. 

However, from the considerations of Sec. Ill B 2 and 
Ref. 14, it seems curious that a direct application of the 
RG equations for the real frequency T-matrix does suc- 
ceed. We want to show here that actually similar (al- 
though less severe) problems do arise. While we believe 
that the frequency behavior given by Eq. (79) is correct 
in the oj — > limit, it turns out that sub- leading terms 
are incorrectly given by the present calculation. To see 
this clearly, we extend our results to the metallic case 
r = 0. When u; ^ T, the integral in (70) behaves as 
I{lj/T) ~ ln(w/T) + ao- aiUj/T with ao ~ 0.125 and 
ai ~ 0.10. Callan-Symanzik equations give in this case: 



-T"iio,T,j>0) 



dynamics for the frequency dependent T-matrix directly 
from the RG equations. 

A generic expression for the critical crossover function 
can be further given by noticing that the renormalized 
T-matrix, at the level of the CS equations, reads: 
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(77) 



so that leading perturbative corrections vanish for the 
choice A = 



--f('^/r)|^|y'^ [note that this alternative 
choice for A interpolates correctly between T/ fi and uj/fJ. 
for both limits of small and large w/T]. At criticality, 
one thus has: 



-T"{iu,T,j = r) = j2^o^e2'-^("/^) (78) 



2r 



which corresponds to the top curve shown in Fig. 11. 
We now turn to the crossover regime (j < r) at finite 
temperature. In this case, the complete result at one- 
loop order reads: 



-T"{u;,T,j<r)^j'No 



U;>2r 



[{T*y + |u;|''e-'-^("/T)]2 



(79) 



which corresponds to the formula for the full scaling func- 
tion ^T{i^/T,T/T*) defined in (17). 

This function is plotted in Fig. 11, and we briefly com- 
ment on this result. First, we notice the appearance of 
an asymmetric hump for T{lu,T,j) at w ^ T, while such 
feature is absent for T{uj,T = 0,j) at w ~ T* (see the 
inset on Fig. 10). This is due to the presence of am- 
plitude corrections in the first case, formula (79), which 
are absent in the second case, formula (71). The shape 



[Hu;/Tk)-I{u;/T)]^ 

1 
[ln(r/TKe"«)-Hai|w/r|]: 



(80) 
r(81) 



This would imply both a shift in the Kondo temperature 
and a low-frequency cusp, which are physically wrong. 
This type of artifact is related to the difficulties which 
come from the use of RG equations for real frequency 
quantities, as discussed previously. While they are quite 
severe in the metallic Kondo problem, the fact that the 
finite temperature low-frequency T-matrix in the pseu- 
dogap case generically vanishes as uj^ at low frequency 
makes such corrections practically hard to see. Indeed, 
at the critical point, formula (78) implies an incorrect 
prefactor e^''"", which leads to a negligibly small overall 
amplitude correction to the exact result. 

We thus conclude this section by emphasizing that in 
the perturbative regime, r <C 1, all scaling functions of 
the pseudogap Kondo problem have been computed with 
reasonable accuracy. 



IV. CROSSOVER FUNCTIONS NEAR r = 1/2: 

WEAK-COUPLING REGIME OF THE 

ANDERSON MODEL 

We want here to investigate the vicinity of the quan- 
tum critical point SCR, now on the side where the model 
flows to the strong coupling SSC fixed point, see Fig. 1. 
As detailed in Ref. 9, we have to use the particle-hole 
symmetric Anderson Hamiltonian as the effective model. 
The expansion is performed around the non-interacting 
resonant-level model, with the Coulomb interaction U 
between the localized electrons being the expansion pa- 
rameter - this allows a perturbatively controlled calcu- 
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lation for DOS exponents r close to 1/2. The starting 
point is the action 

S = -- ^ da(wn) MoSgn(w„) |w„|''](i^(w„) 



+ I'dr U (d,d, - 1) (d,d, - \ 



(82) 



where the local d^ fermions are dressed by conduction 
electrons (u;„ above denote fermionic Matsubara frequen- 
cies (2n + l)7r//3). Technically, this expression can be 
obtained by integrating out the electronic environment 
in (3), giving Aq = ■kV'^NqCOS~^ (t")- "^^^ interaction 
term is written in an explicitly particle-hole symmetric 
way. We can define a dimensionless coupling u, related 
to the bare coupling [/ by w = A^ Z~^ Zfj.'^U, where /i is 
the renormalization energy scale, and we have absorbed 
the non- universal number ^o- Further, — e denotes the 
bare scaling dimension of the coupling U, with e = 1 — 2r. 
The renormalization factors Z4 and Z are associated to 
the vertex and field renormalization. The beta function 
to two-loop order was explicitly derived in Ref. 9: 



P{u) = eu - Au^ , A ■ 



3(7r - 21n4) 



(83) 



with the fixed point u^ = e/A. The renormalization fac- 
tors read 



Z4 = 1 



Au^ 
2e 



Z=l, 



(84) 



the latter relation being exact. ^ 

In a manner similar to Sec. Ill we now calculate the 
running coupling constant, which is simply obtained by 
the substitution k, — v?. We find 



u\\) 



e/A 



1 + [(2e - u)/u]\'- 



2e 



(85) 



We can write this coupling in a similar manner as 
Eq. (22), yielding for u < Uc = \/e/A 



u'{X) = 



e/A 



1 



(^) 



-2€ 



(86) 



where the crossover scale from SCR to SSC is given by 

„,2 \ l/(2e) 



T* = fi 



uf. 



(87) 



Note that the "distance" between SCR and SSC is of 
order -^/e (and not e), which makes the present expansion 
numerically less accurate than the one performed from 
the Kondo model. 



A. Impurity susceptibility 

We start by describing the crossover in the static quan- 
tities. It is worthwhile noting that although the criti- 
cal point SCR is identical to the one found using the 



^ Ximpi^ ) 




0.08- 



0.06 



T/D 



FIG. 12; (color online) rximp(r) at r = 0.47 capturing 
the crossover between SCR at high temperature and SSC 
at low temperature. The NRG data (dots) are obtained for 
U/D = 1, V/D — 0.845, while the analytical curves are plot- 
ted according to the perturbative result, respectively the full 
one-loop result (88) for the upper curve, and the strict 0{^) 
result (89) for the lower curve (see text). The extracted 
crossover scale is T* /D ~ 1 x 10~^^ . Note that the agreement 
is as always better at low temperature where the calculation 
is perturbatively controlled. 



weak-coupling expansion within the Kondo model, the 
crossover is now different, since we here expand around 
the strong-coupling fixed point of the Kondo model. Fo- 
cusing on the impurity susceptibility, and using the per- 
turbative expression derived in Ref. 9 together with the 
running coupling (86), we find that 



TXimp(T) 



(l-2-i-'^)2C(l + r)^ 



27r2+2- 



0.19 



1 



1 






u(T/ fi){88) 
:. (89) 



where r ~ 1/2 has been inserted in the prefactor of 
eq. (88) to obtain the last equation. This result is com- 
pared to the numerical data obtained by NRG on Fig. 12. 
We have to note that the present expansion also differs 
on a quantitative level from the one based on the Kondo 
model in the sense that the r-dependent prefactors, such 
as the one in equation (88), vary quite rapidly with r 
(this is not the case for the coefficient in equation (24)). 
Although a strict expansion at order ^/e requires to ex- 
pand the prefactor at order e° as in equation (89), which 
gives the lower curve of figure 12, a much better agree- 
ment in the range r = 0.4 . . . 0.5 can be reached by using 
the full value of the r-dependent prefactor, as seen in the 
upper curve of figure 12. 

As in Sec. Ill A 2, the impurity entropy has tiny varia- 
tions, and we do not study this quantity further. 



B. Local susceptibility 
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Similarly to the calculation done with the Kondo 
model, we write down the CS equation for the local sus- 
ceptibility in the Anderson model. The renormalization 



factor Z^ was calculated in Ref. 9, with the result 



Z. 



2u 
1 

7re 



giving the anomalous exponent 
r\Au) = - 



(90) 



(91) 



Note that 77' differs from r\^ defined previously as the 
expansion is done now around a resonant level and not 
around a free spin limit. These coefficients are related by 
the relation 1 — ry^ = — e — 77' correcting a misprint in 
Ref. 9. The scaling equation for the local susceptibility 
within the Anderson model naturally reads: 

Xxo.{w,T,u) = A^gxWxioc(y,|,w(A)j (92) 



with g^{\) = exp I / 



/3(u') 



du 



(93) 



using the fact that x\o<J\iv^T,u) has bare scaling dimen- 
sion e. The exponential prefactor can be evaluated to: 



ffx(A) = 



(u(A) - y^) (u + v/^) 



(94) 



The perturbative expansion for the local susceptibility 
is given in diagrammatic form in Fig. 13, and simply 
reads at lowest order in u: 



Xlociil^n) 



[XQ{il'n)+UH "Xoiit^nf] 



(95) 



In the above expression, xo is the susceptibility of the 
resonant level, and involves the computation of a simple 
bubble graph. (Again, the renormalization factor Z4 is 
unity to the order needed here.) Focusing first on the 
T = case, we find Xoi'i'^) = —W\'^/{'^^) as e ^ 1. We 
remark that a direct evaluation of the Z^ factor from 
Eq. (95) and the above expression would provide the in- 
correct result Z^ = 1 — u/(TTe). This is because Xo{'^^) 
presents already a pole at bare level, invalidating the es- 
timation of Z^. To do things in the proper way, it is 
better to perform the analytic continuation in Eq. (95) 
and take its imaginary part: 

XU^) = Z-W^{cj) [1 + 2u^,-^Xo{^)] (96) 

where the real and imaginary parts of the bare suscepti- 
bility at T = are given by the expressions 





FIG. 13: Diagrammatic contribution to the local suscepti- 
bility of the Anderson model up to order u, where full lines 
denote the propagator of the hybridized d level, and squares 
the U — (i'^AqU vertex. 



(above e <C 1 has been used to simplify expressions). 

Note that now x'ai^: T = 0) is well defined in the small 
e limit and that an extra factor 2 appears in (96) with re- 
spect to (95), which is crucial to finally obtain the correct 
Z^ = 1 — 2u/{ne) factor. 

Finally, the CS equation with the choice A = uj/fi re- 
sults in: 

x'L{^,T = 0,u)^^g^{Lo/t,)\u\%gniu) (99) 

which cr osses over from to'^ behavior ai uj ^ T* to 
^e-(2/w)^<i/A bgjiavior at lu ^ T*, in agreement with 
the anomalous exponent value'^^ 

''^ IT A 



(100) 



The above complete expression (99) is compared to NRG 
results in Fig. 14. 

Next we turn to the finite-temperature behavior of the 
susceptibility. As opposed to the case of the Kondo model 
where xioc starts only at order p, the non- interacting 



uj/T-" 



10' 



10- 



10" 



(>oe^ oo°Q°Q 



10 



10 



c(-'^=0) 



10^ 



10- ^^-40 




X^(c^,T = 0) = -^cos(W2)^-^ (97) 

Xo(^,7^ = 0) 



sin(^£/2) \lo\' \ojY , . 

— -^— sgn(c.)c.— sgn(c.) 



(98) 



FIG. 14: (color online) Imaginary part of the dynamic local 
susceptibility x'Ld'^i T = 0) at zero temperature for r — 0.47. 
The analytical result is in Eq. (99). 



16 



u = limit already shows a non-trivial result for xioc 
that includes a temperature dependence as well. Com- 
puting the bare bubble at T > 0, we indeed find: 



+ 00 



r^ / ' ™ _ /x 

Xoi^,T) = -^ / dx x^'' t&nh l^- 



■ 


UJ 




— r 


/O) 


\ 


LU 




— - 


- X 


sgn 





-x] - 


— 




T 




\T 




T 



XU^,T) = — y dxx-'-tanh(- 



(101) 



(102) 



This expression recovers the zero-temperature result (98) 
when Lo ^ T, and gives in the opposite limit: 

x;,(co«r) = -—+T'-Xog— (103) 

Tie TT 



■ 


OJ 




— r 


UJ 


— r 




— - 


- X 


— 


hx 






T 






T 





Xo(^ « T) 






(104) 



Applying the CS equation in this regime (for simplicity 
we focus on the critical point u = Uc) gives for the sus- 
ceptibility (choosing A = T/^): 

\ l-i / ^7re 1 n 

(105) 
which is the expected relaxational behavior. In the 
present case, the success of the RG is due to the fact 
that the bare susceptibility already presents the correct 
small Lo/T limit. 




FIG. 15: Self-energy of the Anderson model at order u^ 



by the self-energy, displayed in Fig. 15, which at zero 
temperature amounts to (for r ~ 1/2): 



E<i(r) = U^Gdoir)' 



-2e 2 
/i U 



1 cos(7rr/2)r(l - r) 

TT \t\ 

1 1 



n3 



-Ai 



(27r)3/2 



Ii-''sgn(r) 
-sgn(T) 



|3-3r 



(107) 
(108) 

(109) 



The Fourier-transformed self-energy reads S(«cij) = 
-ilj,^'^^{u'^ /'iT)\uj\'^^^'^sgn{uj), so that the propagator is: 



Gj^^iiu) = Gj^^(iuj)-T.d{iuj) 



i\uYsgn{uj) 



W- 


LJ 


2e" 


TT 


M 





(110) 

(111) 



The analytic continuation gives for the pcrturbative 
propagator at order u^ 



G'^{uj,T = 0) = -2V2|t 



1 



1 + U^/tt 



(112) 



We are now in the position to apply the CS equation 
with A — Lu/ fi (there is no multiplicative term here as 
Zj- = 1), resulting in: 



C. T-matrix 

In the Anderson model, the T-matrix is simply propor- 
tional to the propagator Gd of the local d level. In the 
non-interacting limit (u ~ 0) one has the free propagator 
Gdo{iijJn) = [isgTi(uJn)\uJnY]^^ whosc analytic continua- 
tion reads: 

G'doiuJ.T) = sm{nr/2)\u;rsgniLu) c^ ^\u;rsgniu;) 

G'^,{u;,T) = -cos{nr/2)\u;rc,-l=\Lor (106) 

using r ~ 1/2 to simplify expressions. Clearly, the T- 
matrix obeys T(a;) ex ^ at tree level (i.e. at the SSC 
fixed point). Since we have already seen in the study of 
the pseudogap Kondo model that the same behavior ob- 
tains at the quantum critical point, we should expect the 
absence of an anomalous exponent (see Ref. 9 for a proof 
of this statement), and simple amplitude corrections for 
the crossover between SCR and SSC (see also Ref. 21). 

The corrections to the bare propagator brought about 
by the local Coulomb interaction u are simply expressed 



T"{uj,T ^0)oc |t^|-'' 



1 



ttA1 + \lu/T*\ 



-2c 



(113) 

This expression is plotted in Fig. 16, together with rc- 
scaled NRG data showing reasonably good agreement. 

We finally conclude on the finite-temperature behavior 
of the T-matrix. From the study of the Kondo model, 
we know that it should obey T"{u},T,Uc) ~ \uj\'' /T^'' for 
o) <C T at the quantum critical point. However, we start 
at tree level with a |w|~'' dependence with r ~ 1/2, so 
that the use of CS equations on the real-frequency axis 
is bound to fail, as one would need a large (hence non- 
perturbative) correction of the exponent. At this stage it 
is unclear whether the strategy of using the GS equation 
directly in imaginary time may be useful here. 



V. CONCLUSION 

In this paper we have undertaken a detailed study of 
the crossovers that occur in the particle-hole symmetric 
pseudogap Kondo and Anderson models, by computing 
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w'^T"{uj,T = Q) 




0.18 



0.16 



Tx 



imp 




T/D 



FIG. 16: (color online) Rescaled T-matrix -uj'^T" (ijj,T = 0), 
Eq. (113), at zero temperature for the same parameters as 
in Fig. 12, showing the crossover on the amplitude of the 
common power law uj~^ . The ringing at low frequency in the 
NRG data comes from numerical errors in the calculated ratio 
related to the broadening procedure. 



FIG. 17: Impurity susceptibility rximp(r) for r = 0.4. The 
upper curve is obtained from the Kondo expansion result 
Eq. (35), while the lower curve originates from the Anderson 
one, Eq. (89). The NRG data (dots) are taken on both sides 
of the quantum critical point, illustrating flow towards the 
local moment and the resonant level fixed points respectively. 



various universal scaling functions associated to the mag- 
netic and electronic response. Our analytical calculations 
are based on the use of the perturbative renormaliza- 
tion group and Callan-Symanzik equations, a powerful 
method that was not previously pursued in its gener- 
ahty for quantum phase transitions. The comparison 
of our results to the available numerical data obtained 
by NRG simulations was seen to be impressive, validat- 
ing the present approach. Some predictions for the low- 
frequency limit of dynamic quantities at finite tempera- 
ture were also made, a regime which remains challenging 
for the numerics. 

As a final illustration, we have compared in Fig. 17 the 
NRG data to our lowest-order estimates for the impurity 
susceptibility, both for small r and small 1/2 — r, taking 
a common value of r = 0.4 which lies beyond the limit 
where our expansions are quantitatively accurate. This 
nevertheless illustrates in a very nice manner how two 
seemingly different expansions are able to capture both 
sides of a common quantum critical point. 

Finally, we think that the formalism that we devel- 
oped and illustrated on these peculiar quantum impurity 
models will be useful for a wider class of problems show- 
ing quantum criticality, e.g., for impurity models with 
bosonic baths and for bulk quantum magnets. 
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APPENDIX A: FIELD-THEORETIC 
DERIVATION OF THE T-MATRIX 



In this Appendix we derive the expression for the T- 
matrix of the Kondo model within a path integral rep- 
resentation and provide a proof of the formula Zj- — 
Z'jZY" for the associated renormalization factor. Follow- 
ing standard literature in field theory ^^, we can express 
the single-particle matrix elements of the T-matrix using 
reduction formulas to the full Green function. Because 
of the local nature of the interaction, this can be cast in 
the simple form^"*'^^ 



Gc{iuJn,k,k') = G'co(ia;„,fc)[l -I- T(icj„)Gco(iw„, A:')] 

(Al) 
where Gco stands for the bath Green function of the free 
system, whereas Gc represents the interacting system. 
We will start by writing out the generating functional 
for the fully interacting Green function of the electrons 
in the Kondo problem^^ 



Zilla,Va) = / ^ 



D [fja 



(A2) 
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where the full action of the problem reads: 
S = -^y^ / d'^kca{iujn.,k)G^Q{iuj„,k)ccr{iuJn,k) 



- -^ / d''/c[rj„c<,(iw„,fc)+h.c.] 



(A3) 



Using this generating functional we obtain the fully in- 
teracting Green function by derivation 



(52lnZ(^,77) 



Sri„{iujn,k)6ri„ii^n,k') 



v^=n„=o 



= Gc{iLOn,k,k'). 



(A4) 



performing the shift c^ -^ Ccr + [G'co?7o-] (this is understood 
as a product in w„ space or a convolution in r space), we 
arrive at: 

S = --'^ d'^kCa{iU}n,k)G'^Q{iljJn,k)Ca{i0Jn,k) 

iuJn 

iuj 

iuJri 

+ j|dr7,^/,,.c„(0)^[ry^(0)Geo] 

+ J jdrjJ-^U, ■ [^jO)G,o]^[r/MO)Gco] 

(A5) 



Taking the derivative of this expression w.r.t. the source 
terms 77 and 77, we find the following expression for the 
fully interacting single-particle Green function: 



Gc{iuj,k,k') = Gca + J^GcQ{T„{iuon)Tl{iun))Gco 

(A6) 



JGcoCy . fl-^fl3'^<r<T)GcO- 



al3 



where we have introduced the composite fermionic op- 
erator To- = Y.a'afj-fl^fp^'yy'^'y'i^)- Comparing this 
expression with (Al), we can identify the T-matrix as 



T{iuJn) 



j2(T,(zc.„)T^(za;„)) + J(^ /t !:^/^t;.) 



a/3 



r = 





(A7) 



+ 



FIG. 18: T-matrix in diagrammatic form: The double solid 
line is the full impurity local propagator of the Abrikosov 
fermions f^, the double dashed line represents the full lo- 
cal propagator of the conduction electrons Ca, and the filled 
square denotes the full vertex function. The second graph 
vanishes in absence of a magnetic field. 



which is illustrated in Fig. 18 in diagrammatic language. 
The second graph, which corresponds to the second term 
in equation Eq. A7, vanishes in the absence of a magnetic 
field. It is now only a matter of identification to note that 
expression for Zr is given by^^ Zr = ZJ^Z^-. 
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